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Q\ . Abstract 

. We compute the two-loop corrections to the coefficient of the b — > 57 magnetic pen- 

guin present in the limits of heavy top and/or heavy higgs. This kind of corrections 
affects in a significant way the observables measured at LEP. On the contrary we 
find that, due to a numerical accident, the correction to BK(B —>■ X s j) is negligible 
(below the 1% level for any possible value of the higgs mass) when the leading order 

CN ' result is expressed in the usual way in terms of the semileptonic BR(I? — > X c ev). 

> ' 

; 1 Introduction 

I The b — > S7 decay seems the more promising observable in B, K and D physics to search for a supersymmetric 
effect. Some observables (like the K — > ttvv decays and various tree level processes) are predictable with 
remarkable accuracy within the SM jjj, but are expected to receive too small supersymmetric corrections. In 
other cases (like in B mixing, or in charm-less B decays) the supersymmetric corrections can be large enough to 
affect the experimental result, but remain masked by too large hadronic uncertainties in the SM prediction 
Qh. Even if the minimal supersymmetric SM (MSSM) Lagrangian contains the new sources of flavour and CP 
violation that can be left at low energy by unification-scale physics , it does not seem possible to disentangle 
the new supersymmetric effects in B-physics from the SM background More optimistic conclusions can 
only be obtained assigning ad hoc values to the ~ 100 parameters of the MSSM Q. or if i?-parity is broken in 
appropriate way. Various supersymmetric effects become larger if a stop state is light (this possibility allows 
baryogenesis at the electroweak scale , but requires an unnaturally large fine tuning g ) or if tan (3 is large 
(this possibility can be realized naturally with some appropriate hierarchy between the parameters of the higgs 
potential). In both cases also the supersymmetric correction to b — > sj is enhanced. 

In all cases, it is interesting to study the b — ► 57 magnetic penguin and it is reasonable to expect that its Wil- 
son coefficient, usually named C7, receives a relatively large and hopefully detectable 0(10%) supersymmetric 
correction. Consequently, it is useful to perform accurate multi-loop computations of this observable. 

In order to obtain a sufficiently accurate SM prediction, the next-to-leading (NLO) QCD corrections, of 
relative order ff§ (a*) /47t, have been computed for all relevant values of the MS scale p: for p near the electroweak 
scale flw M, for p near the B scale pb [|l0|, and for p between pw an d p^ [JTJJ. All these three stages give 
numerically relevant QCD corrections of order 20% to the decay rate. Also the supersymmetric contributions to 
C7 are affected by relevant QCD corrections. However, only the QCD corrections to the charged higgs mediated 
contribution are known |T^| . It is difficult to present in a compact and accurate form the SUSY-QCD correction 
to the remaining relevant supersymmetric contribution, mediated by charginos and squarks. 

Electroweak corrections, of relative order 3|/(4 7r ) 2 an d ln(/2^//2^) • e 2 /(47r) 2 , have partly been computed 
in |ll| and found to be relevant. 

In this paper we compute, for arbitrary value of the Higgs mass, the heavy top corrections of relative order 
<7 2 /47r where 174 is the top quark Yukawa coupling in the SM. We also compute the corrections of relative 
order A 2 /Air where A is the quartic higgs coupling in the SM. More precisely the couplings gt and A induce 
corrections that depend on the top, higgs and W boson squared masses (m 2 , m^, m 2 ^), and increase when mt 
or nth become large, as summarised in table [l]. We compute these potentially relevant corrections in the limit 
mw — ► 0. In the same approximation, the 'heavy top corrections' and/or the 'heavy higgs corrections' to 
the precision electroweak observables measured at LEP have been computed in (m) and/or p5|. The heavy 
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top corrections to LEP observables are so relevant that also the sub-leading terms, suppressed by a power of 
m w/ rn t' have recently been computed pq| . From these results it is possible to derive the heavy top corrections 
to various B and K decays jl7| generated by 'electric' effective operators. On the contrary, a new computation 
is necessary to obtain the corrections to the b — > s"f magnetic penguin. 

This computation is outlined in section 2 (various details are confined to the appendices) and is more 
cumbersome than in the other mentioned cases. It is not possible to relate the magnetic penguin to a simpler 
vertex, as in jlj]. As a consequence more than 50 two loop diagrams, shown in fig. |^, need to be evaluated. 
Rcnormalization is a complicated task: we cannot set the 6-quark mass nib to zero, we have to deal with CKM 
mixing, and we have to connect b — > S7 with the measurable B — ► X s j. A particularly appropriate way of 
doing the renormalization is described in sections 2.3 and 2.4. The final result is discussed in section 3. We 
find that the percentage correction to BR(i? — > X s ~/)/ BR(B — > X c eu) is negligible — at the 1% level — for 
any interesting value of the higgs mass. This is due to a numerical accident; the single contributions would 
separately give few % corrections. The presence of accidental cancellations suggests that the approximation 
mw — ■* could be not a good one. In any case it appears unlikely that the neglected corrections, suppressed 
by powers of m^/mf, be sufficiently large to make the heavy top and/or heavy higgs effects relevant. 



2 Computation 

In this section we outline the computation. For simplicity we will refer to the heavy top and/or heavy higgs 
corrections that we want to compute as l O(gt, A) 2 corrections', even if this is not exact. The relevant Lagrangian 
is given in appendix A, where the various parameters are defined. 



2.1 One loop result 

We begin with recalling the structure of the one-loop result. The contribution of the first generation is negligible, 
so that the unitarity constraint on the CKM matrix becomes V c bV* s + VtbV t * s — 0. It is also unnecessary (even it 
it would be immediate) to compute the 0{gt, A) 2 corrections to terms suppressed by the ratio m s /irib between 
the masses of the strange and the bottom quarks. The effective Hamiltonian for the b — ► sj decay is 

r 8 

n ee = -^v tb v t * s J2 c i°i C 1 ) 

V 1 i=l 

where Gp is the Fermi constant, Cj are the Wilson coefficients and, for our purposes, the relevant operators are 
2 = 4(s Vlc 1 )^^), 7 = 4 7 /^m b (s l7 ^ R A)iV, = Ijj^mbisi^TfM^G" 

[iTTJ [iTT ) 

where i,j are colour indexes. At leading order, and before including QCD corrections, C2 = 1. The leading 
contribution to C-j is given, in the the renormalizable Fey nman- Fujikawa gauge |Q, by the one- loop graphs 
(named "W" , "I" and "E" ) shown in fig. [p 

C 7 = [P w {x t ) - P w (x c )] + [ - ^Pi(x t ) + Y P E(xt)] = \xtP E {xt) « -0.2 (2) 

where Xt = mf/m^, x c = rrvjrrv^ w and the three penguin loop functions P(x) are given in appendix C. In 
the limit m t ^> m\y (in which we will compute the 0{gt, A) 2 corrections) Pw(x t — > 00) = and the one-loop 
result (H) simplifies to 

C?(m t » m w ) = [0 - g] + [1 - 1] = -1. (3) 



2.2 Two loop 0(g t ,X) 2 graphs 

We have seen that the one-loop b — > 57 decay, in the limit of heavy top, does not receive any contribution from 
the graphs "W" of fig. |l| with a VF-boson and a top quark in the loop. It is easy to understand this fact with a 
simple dimensional analysis and a look at the vertices, explicitly written in fig. |l|. For the same reason 01 also 
at two loops the heavy top and/or heavy higgs corrections can be computed in the gauge-less limit. In the same 
way the contribution mediated by the W boson and the charm quark does not receive heavy top and/or heavy 
higgs corrections. At the light of these considerations, all the two-loops graphs that give 0(g t ,X) 2 corrections 
are shown in fig. || (the 12 graphs named 'B', '/i' and 1 xl vanish in the limit of zero higgs mass). We have not 
plotted other 10 graphs of tadpole type (we will discuss them in the following). 
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Figure 1: The one-loop b — ► 57 graphs. The symbol ® denotes a m b factor. All possible attachments of the 
photon have been shown. 



All the graphs contain a (bottom quark 6)-(top quark i)-(unphysical charged higgs x) vertex. The SM 
Lagrangian ( |A.3| ) contains two different vertices of this kind: 

9bb R t L x~ and - g t b L t R \~ ■ 

With the exception of the 1 Kl\ 'r/x,' and l \h graphs (marked with a <g> in fig. |^) both vertices contribute to the 
b R - * sl7 magnetic penguin, and give two different loop structures. The difference is clear from the one- loop 
graphs of fig. [j], where we have separately plotted the two contributions, named "I" and "E". The "E"-type 
graphs are usually more cumbersome to compute, give smaller contributions and generate a different b — > S7 
operator, equivalent to O7 only on-shell. 

In total it is necessary to compute 64 two loop graphs. This computation can be done in few minutes using a 
Mathematica Jl9| code. Since we never need to define traces like Trjsj^ • • • 7 P , nor the completely antisymmetric 
tensor, we can employ naive dimensional regularization (i.e. anticommuting 75) with MS rcnormalization scale \x. 

2.3 Renormalization procedure 

Unrenormalized tree level quantities will be denoted with a supscript °. The renormalization factor of a field <j> 
(cf> = {tL, tn, bL, bn, X, ■ ■ ■}) is named Z$, and <fp — <j)j ' ^J~Z$. The renormalization factor of a parameter p (p = 
{m t , m b , v, V ts , ■ ■ •}) is named Z p , and Z p p° is p° plus the quantum corrections to it. All the renormalization 
Z factors arc defined in a more precise way in appendix B. 

Renormalization would be a very cumbersome task if done adding counterterm diagrams to cancel the one- 
loop subdivergences. It is more convenient to proceed in a different way. Renormalize amounts to express the 
quantum-corrected b — > sj amplitude in terms of quantum-corrected physical quantities: the bare parameters 
must be substituted with measurable quantities, and the b, s, 7 fields have to be re-normalized. In the limit 
n%t raw the one loop b — > 57 amplitude depends on the following bare parameters: g® (or mf), g® (or m°), 
v° (or Gp, that appears as the overall dimensional factor in ([fl)), the electric charge e (that do not get any 
0{gt, A) 2 correction, because it must be the same for any quark or lepton) and the CKM mixing angle between 
the two heavy generations. We define the rcnormalizcd parameters in terms of precisely measured quantities: 



the pole top mass, m t , given by 



m t = m° t ^= K= 5t V) (4) 

V Z tL Z t R 



In the limit m t ^> mw the one loop contribution to C7 depends on m t only through terms that vanish as 
e — * 0, so that only the divergent parts of the Z give a contribution. It is thus irrelevant to distinguish 
between the top pole mass and top running mass parameter. 

the pole bottom mass, m&, given by 

m b =m° b Zm » «=9° b v°) (5) 

Since the one-loop O(g^) correction to the b propagator is mediated by a heavy top, these is no difference 
between pole (on-shell) and zero-momentum bottom mass. 

the Fermi constant measured in /j, decay, Gp. Including the 0(gt,X) 2 corrections, Gf (or the W mass 
term at zero momentum, mw) is linked to the higgs vacuum expectation value (vev) v by 

V2 2m 2 w v*Z x [0) 
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Table 1: Asymptotic dependence of the heavy top and of the heavy higgs corrections to a typical observable 
measured ad LEP, and to the coefficient of the b — > sj magnetic penguin. 
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Figure 2: The two loop-diagrams that give 0(gt,X) 2 corrections. The thick (thin) continuous lines represent the 
top (b and s) quark. The dotted (dashed) thin lines represent the charged (neutral) unphysical higgses x (v)- 
The thick dashed line represent the SM higgs boson h. The thin wavy lines represent all possible attachments 
of the photon. In the graphs marked with a * the correction on the x propagator has to be subtracted at zero x 
momentum. Tadpole diagrams are not shown. 
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where Z x is precisely defined in appendix B.3. The relation (J^), valid in the gauge- less limit, can be 
obtained juj considering the boson vectors as external currents and noticing that the Ward identities of 
the SU(2)l global symmetry imply that the effective Lagrangian contains the following term 



•£ e ff = Z x 



d^x~^V2W^ + ■ ■ ■ = m 2 w \W\' 



The higgs vev v must be determined from the minimisation condition of the effective potential. In the 
gauge-less limit this is equivalent to require that the unphysical higgses \ an d V be masslcss Goldstone 
bosons: 

2 2 

m x = /i 2 — 2\v 2 + quantum corrections of order ^ ^* and smaller = (7) 

At this point we encounter a technical problem: there are corrections to v 2 of order gf v 2 . When vq is 
expressed in terms of v, the term in the one loop coefficient of the b — > s-f penguin suppressed by one 
power of x^ 1 — m 2 v /m 2 , Cj = —1/3 + 3(lnxt)/4xt + • • •, would give a heavy top correction, infrared (IR) 
divergent in the limit mw —> 0, that cancels a similar IR divergence present in the two loop graphs. To 
avoid these infrared problems, it is more convenient to cancel the quantum corrections to the Goldstone 
masses in eq. (0) adding appropriate counterterms in order to obtain v = v°. This procedure requires to 
include the new b — > sj Feynman graphs that contain these counterterms. This is equivalent to compute 
the two-loop graphs named "F" and "B" and marked with an * in fig-p] with the prescription that the 
correction to the \ propagator has to be subtracted at zero \ momentum^ When this is done, all infrared 
divergences disappear and the graphs can be safely computed at mw — 00- 

Finally we have to renormalize the CKM mixing angles. Since only the mixing between 2 nd and 3 rd 
generation is relevant, the mixing is described by only one parameter. Furthermore we may set V t b = 1, 
Vts = —Vcb = S23 <C 1: we can neglect the 0(gt,X) 2 corrections to terms suppressed by powers of the 
Wolfenstein parameter A 2 ^ 0.05 p0|. At this point, we remember that normalising the b — ► 57 decay 



rate to the semi-lcptonic b —> cev one allows to eliminate a large hadronic uncertainty coming from a m\ 
phase-space factor. For this reason, we will extract the renormalized value of S23 from BR(B — > X c ev). 

We now show how is it possible to accomplish the last point. It is very convenient to renormalize a non standard 
form of the bare Lagrangian that is respected by the quantum corrections that we are considering. Keeping only 
the bare interactions and the bare fields necessary to describe flavour mixing, we write the SU(2)^-symmetric 
bare Lagrangian as 

£0 = £ fipq° + 9° t t° R Q° 3 H*° + g° b b° R (Q° - s%Q° 2 )H° (8) 

where q = Q%, t R , b R }. The down components of the quark doublets, Q?, are not mass eigenstates. Had we 
started from the usual flavour basis of mass eigenstates, complicated flavour redefinitions would be necessary to 
obtain a unitary renormalized CKM matrix. The tree level Lagrangian and the divergent part of the quantum 
corrections have a SU(2)l global symmetry. This symmetry is useful to understand some cancellations (for 
example why Z mt is not divergent) and to fix some dangerous signs. It is easy to verify that the 1PI effective 
Lagrangian, corrected only by the divergent part of the quantum corrections is 

^effective = J2 Z «? i] Pl + Z 9t 9t t R Q°sH*° + g° b b° R (Z gb Q° 3 - s° 23 Q° 2 )H 
1 

and has the same form as the bare Lagrangian (||) (its last term does not get 0(g 2 ) corrections). The Z factors 
are 

Z Q 3 = 1 + 2 (4^)2 e' ZtR = Zgb = 1 + (4^)2 i' Z Q-2= z g t = z bn = l - 

We are now ready to renormalize also the finite part of the quantum corrections. Performing, in the effective 
Lagrangian, the following redefinitions of the fields 





sl - S23 bL 



and couplings 



1 For this reason we have not plotted the graphs with a tadpole correction to the \ propagator: they are cancelled by this 
procedure. The tadpole graph in which the external photon is attached to the tadpole loop also gives a zero contribution to the 
b — » S7 magnetic penguin. 

2 We have checked that the more cumbersome procedure of keeping myv as an infrared regulator gives the same final result. 
When doing the computation at = 0, it would be useful to know if the 1/e poles are due to infrared (IR) or to ultraviolet 

(UV) divergences. However some of the graphs in fig. ^contain terms like e/e uv e; r , so that we have not been able to distinguish 
IR from UV divergences using dimensional regularization only. 
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vg t =m t , vg b =m b — , s 23 = — ==s 23 (9b) 

the renormalizcd parameters mf , mi, and S23 coincide with the physical quantities previously chosen (the meaning 
of the Z factors should be clear from the notation; they are precisely defined and computed in appendix B). 
Notice that, only at this stage, we have performed a unitary flavour rotation in the down components of the 
quark doublets Q, so that the renormalizcd quark fields are mass eigenstates. The V c b CKM mixing angle that 
appears at the gauge cifflbi, vertex is S23 (the gauge interaction of Q% does not get 0(gt,X) 2 corrections), so 
that the link between B — ► X s j and the semi-lcptonic decay B — > X c ev decay can now be easily extracted. 



2.4 Renormalization of C 7 

We are now ready to express the one-loop b — > S7 magnetic penguin in terms of renormalizcd parameters and 
fields. The magnetic penguin obtained from the (unusual) tree level Lagrangian (||) is 

where C?[a, 6] = -pfp- (07^ 7^6)^1,. Notice that no explicit mixing angle appears at this stage. The two 
contributions arise respectively from the one-loop graphs named "I" and "E" in fig. |l| and give rise to two 
different operators, that are equivalent only on-shell. Since the graph "E" is not one-particle irreducible, we 
must take this fact into account in the renormalization. Operating the redefinitions (^) and (||), and using the 
renormalizcd Dirac equation idb = m^b to set the 6-quark on-shell, we obtain 

<T„,,. - - lit} + <"> 

The final formula for the coefficient C7 of the b — > S7 magnetic penguin, that includes all renormalizcd 0(<?t, A) 2 
corrections, is thus 



23 5 Z x 1 Z x g 2 If /two- loop diagrams\ 
7 ~~ ~36 + 12^~~~ ~ 9 ~ (4tt) 2 24 + \plotted in figure 2 J ' 



(12) 



We have added the charm contribution, —23/36, that does not get 0(gt,X) 2 corrections. In conclusion the 
heavy top and/or heavy higgs correction to C7, computed in the limit raw — * 0, is 



with 

cf = 



C 7 (ra 2 » m 2 w ) = -1 + j^cf* 



-16 + 39r - fir 2 - 26r 3 -16 + 2r - 36r 2 + 74r 4 - 45r 5 + 2r e 



rr 2 



144r 864r 2 

8 - r - 6r 2 - 52r 3 + 85r 4 - 33r 5 + 2r 6 . , 80 + 68r - 262r 2 + 134r 3 - 25r 4 + 2r 5 ,r\ 

Li 2 (l - r) <&(-) 

72r 2 V y 288r V 

8 - 17r - 2r 2 - Ur 3 , 74 - 45r + 2r 2 



lnr - — — — r 2 In 2 r w 0.14 - 0.0024r - 0.046 lnr 



72r 288 

where r = ?n 2 /m 2 , the functions Li 2 and $ are defined in appendix]^, and the approximation holds for any 
reasonable value of rat G [60, 1000] GeV. Since, instead of using MS subtractions, we have expressed the result 
in terms of physical quantities, no dependence on the MS scale p. is left. In the limits mu *C ra t and 3> ra t 
C-j reduces to 

2 7T 2 1 7T 2 1 2 1Q 17 7T 2 

C?* fr 0) = AT C (— --)-(—- — ), C* 7 9t (r -> 00) = -— lnr + — + — 
7V ' cV 36 9 ; v 36 16 ; ' 7V ' 144 432 24 

where N c = 3 is the number of colours. In table [l] the asymptotic dependences of these corrections in the limits 
of heavy top and of heavy higgs are compared to the corresponding limits of the corrections to the observables 
measured at LEP. 



6 



?- 

X 

T 

a! 
cq 

o 

a 

_o 

o 

o 
aj 
CO 

w 

a 

oj 
o 

o 

Oh 




Higgs mass in GeV 



1000 



Figure 3: TTie percentage correction to BR(_B — 
Higgs mass in the limit of heavy top ad/ or higgs. 



X s j) for mi (pole) = (175 ± 5) GeV and as function of the 



2.5 Branching ratio 

With our renormalization scheme the B — > X s j branching ratio can be linked to the b 
the usual way (see jllj and references therein), 
corrections 



that we now briefly recall. Neglecting, 



BR(B -> X s7 ) 



where BR(5 



r(6 — > cez/) 

X c ei^) = 0.105 ± 0.05 has been measured |2l) , 



BR(5 -> X c e^) 



-> S7 decay width in 
for simplicity, QCD 



(13a) 



T(b^ si) 
T(b — > cei/) 



32tt 3 
32^ 



rogl 



S23C23) 2 ^|C 7 (^) 



; TO h( s 2 3 ) 



6 : 



(13b) 
(13c) 



and <?(z) is a phase space factor |TT|. The 0(g t ,A) 2 corrections to the the semileptonic width are entirely 
contained in our definitions of the renormalizcd parameters. As previously said the complete CKM mixing can 
be reinserted neglecting the 0(g t , A) 2 corrections to the terms suppressed by powers of A 2 ^: (S23C23) 2 / (S23) 2 — * 
\Vt b V ts /V cb \ 2 = 1 + A 2 v (2p w - 1) + 0(A W ), (A w ~ 0.22 and p w are Wolfenstein parameters |o)). 

These expressions receive important QCD corrections (pcrturbative and non perturbative) , that can be 
added in the usual way. Infact, the 0(<?t, A) 2 corrections that we are considering, only affect the values of the 
Wilson coefficients at the electroweak scale, Ci(flw), i = !->■■■,&■ We have computed the correction to C 7 . 
Other relevant 0(gt, A) 2 corrections can be present in C2 (the coefficient of the (s 1 '7 '^T i lc*)(c j '7 '^T 'x,W) operator) 
and in Cs (the coefficient of the b — » sg chromo-magnctic penguin). In a formal expansion in powers of a.3, the 
corrections to C2 and Cs enter at the same order as the correction to C7. 

X s 7 decay, in view of the small 



• However, the correction to Cs cannot have a significant effect on the B 
mixing coefficient between C7 and Cs: 



C 7 (flb) « -0.155 + 0.65 C 7 (p. w ) + 0.085 C s (fl w ), 

where ~ mj, and p.\y ~ Tfiw ^ n an Y case, it is immediate to obtain the renormalized correction to C's 
from a subset of the two-loop graphs of fig. ||. In the heavy top limit Cs has a weak dependence on the higgs 
mass rah and can be approximated with its value at = 0, Cs(p-w) = ~ l/8 + .9 2 /(47r) 2 (7r 2 /12 — 13/32). 



• The 0(gt, A) 2 correction to C2 vanishes in our renormalization scheme. 
We remember that the coefficient C7 can also be extracted from the spectrum of the B 



X s £+£~ decays | 
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3 Numerical result 



Due to a numerical accident, the potentially relevant heavy top correction to the b — > magnetic penguin 
turns out to be negligible, if the leading order result is expressed in terms of the renormalizcd parameters Gp 
and BR(_B — > X c ev). The percentage correction to BK(B — > X s ^f) is plotted in fig. || as function of the higgs 
mass, rrih, and for m t = (175 ± 5) GeV, where mt is the pole top mass. In the plot we have included the QCD 
corrections at NLO order and all the heavy top and/or heavy higgs corrections (also the small ones, not fully 
given in the text, due to the mixing with the chromo-magnetic penguin). 

We see that the correction is never larger than 1% for any possibly interesting value of the higgs mass. 
We remember that small values of the Higgs mass are preferred. The electroweak precision measurements give 
rrih < 250 GeV at 90% C.L. [p2| . Moreover, the minimal supersymmetric extension of the SM (MSSM) predicts 
a light higgs: at tree level it requires that rrih < mz- Adding the one-loop corrections |23j and barring an 
unnatural fine-tuning larger than 20, the MSSM bound on the higgs mass becomes m/j < 120 GeV [pH , 

As said, the smallness of the correction is due to a numerical accident. Before having done this computation a 
relevant heavy top correction, at the 5% level or even larger^], could not be excluded. However, due to accidental 
cancellations between different contributions these corrections are smaller than 1%. 



4 Conclusions 

We have computed the heavy top and the heavy higgs corrections to the Wilson coefficient of the b — > 
magnetic penguin in the limit mw — ► 0. These potentially relevant corrections turn out to be very small. Since 
there are accidental cancellations between different contributions, it is possible that the limit — * be not 
a good approximation. However it appears unlikely that the neglected corrections, suppressed by powers of 
rriyy I m 2 , be sufficiently large to make the heavy top effects relevant^ For these reasons it seems safe to conclude 
that these kind of corrections can be neglected. 

Acknowledgements I thank R. Barbieri for several illuminating discussions and M. Ciuchini, G. Degrassi 
and G. Giudice, that collaborated in an early stage of this work. 



A The relevant Lagrangian 



As explained in the text we can work in the gauge- less limit. In a flavour basis where the Yukawa matrices g n 
and gd arc diagonal, the relevant Lagrangian is 



C = qidq + \dH\ 2 + n 2 \H\ 2 - A|ff| 4 - (u R9u QH* + d R g d V^QH + h.c.) 



(A.l) 



where V is the CKM matrix, d = {d, s,b}, u = {u, c, i}, q = {u, d}. The contractions of the SU(2) l doublets Q 
and H, expressed in terms of their up (|) and down (J.) components, are 

urQH* = u R (Q r H^ + QiHD and d n QH = d fl (Q x i? T - Qt^x) 

so that the mass terms of the quarks have conventional sign. Decomposing the fields into mass eigenstates (we 
denote the unphysical higgses with greek characters) 



H 



v + {h + irj)/y/2 

X- 



(u L 



one obtains 
where 

£frcc 



C = £f rco + £yuk - V cu bic - V quar tic + V2v(fl 2 ~ 2Xv 2 )h 

(Oh) 2 (Or,) 2 



(A.2) 
(A.3) 



q(i& - m q )q+ \dx\ + 



(6Xv 2 - f)% + (v 2 



2Xv 2 )(^ + \ X \ 2 ) 



■Cyuk — — 9t 



V2 



9b 



b^—^-b-(x- V tb b R t L +h.c.) 
v2 



3 Even if gzil^w) ~ 1-2 is not significantly larger than gt ~ 1.0, the well known NLO QCD effects of relative order g|/(47r) 2 give 
significantly larger corrections (at the 20% level) because they have 'more colour' and 'more spin' circulating in the single Feynman 
graphs. 

4 The ones enhanced by a N c factor have recently been computed in ]13| . Including QCD corrections they give another — 1% 
correction to BR(B -» X S ^)/BR(B -+ X c ei>). 
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V cu b 1( 



V, 



quartic 



V2Xv(h 3 + hr] 2 + 2h X -X+) 

*X-X+(h 2 + rf + X-X-) + \{h 2 + ?? 2 ) 2 . 



B Renormalization constants 

In this appendix we give the Z constants needed for renormalization. 



B.l Bottom quark 

We write the inverse 1-loop propagator of a b quark with momentum pt as p b [Zb L Vi, + Zb R Vu\ — Z mb rrib where 
Pl,r — (1 T 7s)/2- Since there are no infrared divergences even in the limit mw = 0, the constants Z depend 
on pi, only trough irrelevant terms suppressed by powers of p 2 /m 2 . Their values are: Zb R — 1, 



Zb, 



9t 



(4tt) 2 



-(- +ln-^) + - 

m 2 4 



2 y e 



and 



Z„ 



9t 



(4tt) 2 



- + In — 7) 

e mi 



B.2 Top quark 

We write the inverse 1-loop propagator of a t quark with momentum p t as p t \Z tl Vh + Z tR Vii) — Z mt m t . We 
need only the divergent parts of the Z constants. These divergent parts do not depend on p t and can be more 
easily computed in the limit of unbroken SU(2)x, global symmetry. We find: 

„2 rn i _s 

Z mt = 1 



9t 
(4tt) 2 



Zt R = l + 



at 

(4tt) 2 



Ztr, = 1 



9? 
(4tt) 2 



1 

2s 



where • • ■ denotes the finite parts that we do not need. 



B.3 Higgs 

We write the one-particle irreducible (1PI) propagator of a unphysical charged higgs x with momentum p x as 
i/[q 2 Z x (p x ) — m^Zm^p 2 ,)}. As explained in the text we enforce the vanishing of Goldstone masses, m 2 x = 0, 
inserting appropriate counterterms (equivalent to the non 1PI tadpole graphs that shift the vcv v to the quantum 
corrected minimum). We only need the constant Z x at zero momentum. At one loop0 



1 



1 



(4tt) 2 



(B.4) 



with A 



-2tV4. 



C Special functions 

Our final result contains the bi- logarithmic function ~L\2(z) defined as Li2(x) = — J Q ln(l — xt)dt/t and the 
function 

V I 

4 Im Li 2 exp (2i arcsin y/z) for < z < 1 



EE i 



^-ln 2 (4z) + 21n 2 ^VpA 



4Li 2 i^P^ 



for z > 1 



The one loop penguin functions used in (|2|) arc 

-23 + 67 a; -50 a; 2 2 - 7x + 6x 2 



Pw{x) 
Pj(x) 
Pe{x) 



36(x- l) 3 
3 — 5 x 3a; — 2 



6{x - l) 2 3(x - l) 3 



6(x- l) 4 
In .r 



a; In a; 



7-5a; 



3a;- 2 



• In , 



36(a;- l) 3 6(x- l) 4 ' 

5 The corresponding expression given by Barbicri et. al in jl^ contains a misprint (not present in their final result). 



(C.5a) 
(C.5b) 
(C.5c) 
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